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I. INTRODUCTION 

The detection of spontaneous Hawking radiation (HR) 
remains a major challenge in modern physics. Originally 
predicted for cosmological black holes [lj, it was soon 
noticed [1 1 that the emission of Hawking radiation 
(HR) is a kinematic effect that can be observed in an 
ordinary laboratory. For a quantum fluid traversing a 
subsonic-supersonic interface (which amounts to a sonic 
event horizon), the spontaneous correlated emission of 
phonons into the subsonic and the supersonic regions 
has been predicted Fil-ITTI|. Sonic event horizons have 
been realized by accelerating a Bose-Einstein condensate 
[3 ]- In a similar setup, the self-amplifying stimulated 
HR resulting from the black-hole (BH) laser effect E3 
has been observed [l3|. However, the emission of spon¬ 
taneous Hawking radiation still remains unobserved. 

In the context of quasi-stationary flow scenarios, it has 
been proposed that a large leaking condensate can pro¬ 
vide a suitable subsonic-supersonic interface where HR 
production could be conveniently observed [3113. Re¬ 
alistic protocols to produce such quasi-stationary flow 
regimes have been investigated recently |T(> . 

It has also been noticed sa that the violation of classi¬ 
cal Cauchy-Schwarz (CS) inequalities EMI by outgoing 
quasiparticles can provide unambiguous evidence of deep 
quantum behavior and ultimately of the existence of a 
spontaneous contribution to HR. The distinction between 
spontaneous and thermal (stimulated) contributions is a 
fundamental requirement in the search for spontaneous 
HR. Proposals based on density-density correlations in 
real space do not meet that requirement [3 Eoj - The 
mere measurement of phonon [2l| or atom [14j| intensity 
spectra would not permit that distinction either. 

An alternative scheme to identify deep quantum be¬ 
havior relies on the detection of entanglement between 
the various outgoing radiation channels from a sonic hori¬ 
zon [3 ES • A recent work has addressed the possibility 
of detecting entanglement through density-density corre¬ 
lations in Fourier space [3] ■ Approaches based on entan¬ 
glement detection have been also proposed in analogous 


contexts such as inflationary cosmology E3> astrophys- 
ical black holes (!26j, general relativistic quantum fields 
[3 | , or other black-hole analogs [3 E3 ■ 

Our present work aims at clearly establishing the the¬ 
oretical relation between CS violation and entanglement 
in the outgoing quasiparticle modes of analog HR, uni¬ 
fying the existent work of Refs. El El, El- We also 
study the ir p otential measurability in specific detection 
schemes [24| . Importantly, by introducing envelope- 
modulated Fourier transforms we explicitly take into ac¬ 
count the role played by the spatial location of the asymp¬ 
totic regions. This allows us to show that the violation of 
only quadratic CS inequalities can be measured. In par¬ 
ticular, we prove that the complete implementation of 
the generalized Peres-Horodecki (GPH) criterion [3 E3 
or the quartic CS violation E2 requires the knowledge of 
parameters that are impossible or very difficult to mea¬ 
sure, at least within the class of detection schemes here 
considered. However, we also show that this does not rep¬ 
resent an important limitation in practice, since, under 
rather broad conditions, all the previous criteria become 
equivalent. 

Although the present work is motivated by the quest 
for the observation of spontaneous HR, the results here 
obtained are also of relevance in neighboring fields such as 
quantum optics [3 El and quantum information physics 
13 l33j , as well as in the broader topic of bosonic con¬ 
densates [3, :35j- Atom flow through a sonic black hole 
may also be viewed as one of the paradigms of atom 
quantum transport through a barrier [3] . In general, a 
subsonic/supersonic interface can behave as a basic cle¬ 
ment that provides novel functionalities in atom circuits 
within the emerging field of atomtronics [3 . (3l. For in¬ 
stance, the black-hole configuration of Ref. jT6[ can be 
used to produce a quasi-stationary supersonic atom cur¬ 
rent. Another recent numerical work |.39l | has shown that 
a BH laser configuration is able to reach a regime of con¬ 
tinuous emission of solitons, providing a hydrodynamic 
analog of an optical laser. 

This paper is arranged as follows. In Sec. El we 
present the physical model that we use in this work. Sec¬ 
tion nn contains a detailed comparison among the exis- 
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tent works on the detection of the spontaneous HR and 
unifies the criteria under certain conditions. In Section 
IIV1 we discuss the possibility of experimentally verifying 
the discussed criteria. Finally, in Section El we present 
numerical results on the detection of the CS violation in 
typical HR setups. 

II. PHYSICAL MODEL 


their corresponding annihilation operator, with the in¬ 
dex i taking values u— ,dl—,d2 — in. For instance, the 
scattering state d2 — in has the asymptotic form 

Zd2 — in,a; (x ^ Oo) — S u d2 (^) Su— out,cj (*^) (4) 

Zd2— in,aj {x ^ Oo) — 5ci2—in,aj (x) 

+ Sdld2 (w) Sdl — OUt,OJ (pt) 

+ Sd2d2 (w) Sd2-out,w {x) • 


We start by considering the stationary flow of a one¬ 
dimensional BEC in a regime known as ID-mean field 
regime HE S3- The condensate is described by the sta¬ 
tionary GP wave function 4'o(a;) = \Jn 0 (x)e ie( - x ' > . We 
can define the local sound speed and flow velocity as 
c(x) = lgn{x)/m\ l t 2 and v{x) = h6'(x)/m, respectively, 
with g the coupling constant and m the mass of the 
atoms. In our convention, v(x) > 0. A sonic BH configu¬ 
ration is that involving two asymptotically homogeneous 
regions such that one is subsonic (the upstream region, 
x -A —oo), with v u < c u , while the other one is supersonic 
(downstream region, x -A oo), with Vd > Cd■ The mag¬ 
nitudes c Ut d, v u ,d represent the asymptotic homogeneous 
values of the sound and flow velocities. Hereafter, we use 
units such that h = m = c u = ks = 1. Throughout this 
work, we follow the notation of Refs. 14, 15l l20| . 

Within the Bogoliubov-de Gennes (BdG) approxima¬ 
tion, the field operator is decomposed as: 

%) = vl/ 0 (x) + <p(x)e ie W = [^(P) + <p(x)}e i0 ^ , (1) 

where we have removed the mean-field condensate phase, 
e l6 ( x ), from the definition of the field fluctuations <p(x) for 
computational purposes. The expression of tp(x) in a BH 
configuration is given by: 

POO 

<p{x) = / dw V'[u/_ in , tJ (*)7/- in (w) 

J 0 J 

+ ^/-in,^(a ; )7/-in( W )] 

/’Wmax 

+ / cL}[u d2 -in,u(x)p d2 _ in (uj) 

Jo 

+ ^d2-in,^(^)7d2-in(w)] . (2) 

The components of the spinor £ 0 , w (x) = 

[u a ,ui(x),v a ,uj{x)] T are solutions to the BdG equa¬ 
tions at a given frequency ui. Here we have adopted 
the convention of working with positive frequencies at 
the expense of introducing the negative-normalization 
channel d2. The index I is summed over the conven¬ 
tional, positive-normalization channels (it and dl); see 
Fig. □ for the mode notation. In the asymptotically 
homogeneous regions, the macroscopic wave function 
has the form 

’Fo(x) = \/n r (x)e lVrX+8r , r = u,d (3) 

where 9 r is a constant phase. Thus, the scattering states 
Zi U (x) are combinations of plane waves (scattering chan¬ 
nels) in the asymptotic regions u and d, while 7 i(w) is 


In this expression, s 0iW (x) is the spinor (plane) wave 
function of the corresponding scattering channel a: 


Sa,u [X) 

U a (u) 
V a {u) _ 

N a 


0 ik a ( uj)x 


\J2-K\w a (w) 
= N t 


Ua{u>) 

V a (u) 


+ [iO - V r k a (w)] 
- [W - V r k a (w)] _ 
= [ 2k a M |W - V r k a (w )|]" 1/2 , 


( 5 ) 


and the corresponding wave vector, k a {u) is given im¬ 
plicitly by the dispersion relation 

n 2 a = (u>-V r k a ) 2 = C 2 r k 2 a + k 4 a /^, ( 6 ) 


where S2 a > 0 is the comoving frequency, defined posi¬ 
tive, and the index r = it, d characterizes the asymptotic 
region corresponding to the scattering channel a. The 
group velocity of the scattering channel a is w a (w) = 
[dk a (w) /dio] and it is included in the definition of the 
scattering channels in order to properly normalize them 
to a given quasiparticle flux. The dispersion relation is 
depicted in Fig. [T|for the subsonic and supersonic asymp¬ 
totic regions. 

The expressions for the other scattering states are sim¬ 
ilar to Eq. ©• They all are characterized by an incoming 
channel carrying unit flux and several scattering ampli¬ 
tudes (elements of the S'-matrix) describing the transition 
to the outgoing modes, which are also assumed to carry 
unit flux. The “out” (advanced) scattering states are re¬ 
lated to the “in” (retarded) scattering states through the 
scattering matrix S for a given value of the frequency ui: 


lu-ont 


C 

u uu 

Sudl 

^ud2 

7dl-out 

= 

£>dlu 

Sdldl 

Sdld2 

_ 'J d2— out . 


Sd2u 

S d 2dl 

S d 2d2 


7u—in 

7dl-in 
_ ^d2 —in _ 


( 7 ) 


Here, as often in the rest of the paper, the frequency 
dependence is understood. 

We note that Eq. © is a Bogoliubov-type relation 
linearly combining destruction and creation operators. 
This is due to the negative normalization of the d2 modes, 
see Fig. [Q Thus, the vacuum of the “in” modes can 
be regarded as a squeezed state in the representation of 
the “out” modes. This non-uniqueness of the vacuum 
is at the root of the Hawking effect itself. In order to 
detect the intrinsic quantum behavior corresponding to 
the squeezed character of the zero-point HR, we focus on 
the correlations between the “out” modes. 





















3 



- 1.0 - 0.5 0.0 0.5 1.0 



k 


k 


FIG. 1: Bogoliubov dispersion relation on the subsonic (left, upstream) and supersonic (right, downstream) sides. The blue 
(red) branches correspond to positive (negative) normalization. Here, d(u) denotes downstream (upstream). The Hawking 
frequency, cu m ax, is the frequency above which no HR can be generated. 


III. CAUCHY-SCHWARZ INEQUALITIES AND 
ENTANGLEMENT 


Several criteria have been proposed in order to distin¬ 
guish the spontaneous from the stimulated (thermal) HR 
or the coherent collective modes. First, it was argued in 
Ref. E2 that the violation of CS inequalities can be 
regarded as an unequivocal signature of the presence of 
zero-point dynamics. Specifically, the following second- 
order correlation function between two given modes i,j 
was considered: 

r ij = {lltfjljli) > 0, (8) 

where 7 j is the annihilation operator of the mode i. The 
expectation value of an operator O is taken over the state 
of the system, described by the density matrix p, (O) = 
Ti(pO). The relation 

Ty < \jTuTjj , (9) 

is a CS-type inequality which is always satisfied in a clas¬ 
sical context. The violation of © is characterized by the 
positiveness of the difference 

o ij = r ij - v'TaTjj > 0. (10) 

We will refer to Eq. m as the quartic CS violation. We 
can also define the corresponding first-order correlation 
functions, 

9ij = (7i7j>. C H = (lilj) , ( 11 ) 

whose associated quadratic CS violations, in analogy to 
Eq. (gOD , are given by: 

^ ij — Cjj 9ii9jj > 0 . (12) 

The inequality \gij\ 2 < gugjj is always satisfied; see Ap¬ 
pendix [A] 


Another possible signature of the quantum character 
of the system is the presence of entanglement between 
two modes i,j. In this work, we follow the convention 
of defining entanglement as the non-separability of the 
state of the system. We say that a two-mode state p is 
separable when it can be decomposed as: 

P =' 52 PnP i n ) ® Pn ] ( 13 ) 


where pb, is a state of the Hilbert subspace spanned 


by the i mode. The use of the GPH criterion [41|, |4 
was proposed in Refs. [22j, 


231 ] to identify the entangle¬ 
ment between two modes within the context of analog 
HR emission. The GPH criterion asserts that the state 
is entangled if Vij < 0, where Vij is the GPH function 
defined in Eq. (1A11D . In the particular case of Gaussian 
states, the reverse implication is also true, i.e., the GPH 
criterion is also a necessary condition for entanglement 

Si. 


Some general remarks about CS violations and entan¬ 
glement are in order. The violation of a quadratic CS 
inequality, Eq. m, is a sufficient condition for the non¬ 
separability of the system [43j. In particular, we show 
in Appendix [~Y1 that the quadratic CS violation is a suf¬ 
ficient condition for the fulfillment of the GPH criterion, 
which in turn is known to imply the presence of entangle¬ 
ment. It can also be shown that the quartic CS violation 
is independent of the entanglement of the system (see 
Appendix 0, so separable states can violate quartic CS 
inequalities and viceversa. 

Now we focus on the specific case of analog HR in a 
BEC. We evaluate the correlation functions of Eqs. © 
and CD for the “out” modes at given w: 


r 7'M = (7i'-out( W )7]-out( W )7i-out(w)7i-out(w)) 

9ij(u) = (bl-out MVoutM) (14) 

Cij ) = (rfi —out (k'O'Tj—out (^)) • 
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Gaussian 

incoherent 

CS4 

CS2 

GPH 

yes 

yes 

• GPH 

• 

<S> GPH 

• NS 

yes 

no 

• indep. NS 

• 

=> GPH 

• NS 

no 

yes 

• indep. NS 

• 

GPH 

• => NS 

no 

no 

• indep. NS 

• 

=> GPH 

• => NS 


TABLE I: Logical relations between the different criteria stud¬ 
ied throughout this work in the various physical cases. The 
three rightmost entries in the upper row indicate the different 
criteria for the identification of quantum behavior: CS2, CS4 
stand for quadratic and quartic CS violations, respectively, 
while GPH stands for generalized Peres-Horodecki criterion; 
NS means non-separability. The two leftmost columns define 
the various physical cases here considered. By “incoherent” 
state, we understand a density matrix which is diagonal in 
the representation of retarded quasiparticle scattering modes, 
each characterized by a single incoming channel; see Eq. m- 
The symbol • stands for the uppermost entry in the corre¬ 
sponding column. The abbreviation “indep.” means that, in 
the three lower cases, the quartic CS violation is independent 
of the non-separability of the system. 


Hereafter, i,j = u,dl,d2. The associated quartic and 
quadratic CS violations are characterized by the positiv¬ 
ity of 0y(w) and Arespectively, as defined in Eqs. 
mm- I n a similar fashion, we study the GPH func¬ 
tion Vij (w) for the i. j “out” modes at the same frequency 
w. For simplicity, in this section we will obviate the Dirac 
delta factors appearing in all equal-frequency correlation 
functions [see Eq. m for a complete expression]. 

One may wonder to what extent the previous criteria 
differ when applied to spontaneous HR. We devote the 
rest of this section to prove that, when the state of the 
system p is Gaussian and incoherent in the “in” inodes, 
the quartic and quadratic CS violations become equiva¬ 
lent to the GPH criterion. Specifically, the requirement 
of incoherent incoming modes can be expressed as: 

(ffi— in (^)' T i 7 — in (^ )} = 0 

( 7 !—= ni(u)5ij8(u - u') ^0.(15) 

These correlation functions fully characterize the state 
provided it is Gaussian. It is shown in Appendix [Bl that, 
for that class of states, we only need 7 parameters for the 
computation of the CS violations and the GPH function. 
The case of Gaussian states is important because, within 
the BdG approximation, the dynamics is Gaussian. 

For convenience, we define the complex vector 


OCi{uj) 


Siu (^) \JWu (A) 
Sidi{u) y/ndi{u) 

Sid2 (w) ^/n d2 (w) + 1 


(16) 


From the previous definitions, it is easy to show that the 
only non-zero quadratic correlations for “out” modes are: 


gw = u\ ■ ar 

gd2d2 = \dtd2\ 2 — 1 
Cld2 = Ol\ 2 ■ CXl , (17) 


where the index / stands for a normal, positive- 
normalization (u or dl) mode. First, we study the cor¬ 
relation between a normal and an anomalous (negative- 
normalization) mode. The case I = u corresponds to 
the proper Hawking effect [hi], H3] and the case I = dl 
corresponds to the bosonic equivalent of the Andreev re¬ 
flection Q. 

As we are working with Gaussian states, we can apply 
Wick’s theorem to compute the quartic correlations as a 
function of the quadratic correlations 

Tn = 2gj I = 2\a I \ i 
Td2d2 = ^9d2d2 = 2 (|o:rf21 2 - l) 2 
Tld2 = \cid2 \ 2 + giigd2d2 

= \(Xd2' a l\ 2 + \ a l\ 2 (\ a d2\ 2 -!)■ (18) 

Therefore, the condition for the quartic CS violation 
0/d2(w) > 0 reduces to the simpler quadratic CS vio¬ 
lation 


Q/d2(w) = Ald2 (w) > 0 (19) 

which, using Eqs. OH) and m , can be rewritten as: 

Wd2 ■ a A 2 > M 2 (|ad2| 2 - 1) (20) 

The —1 within the second bracket results from the 
anomalous character of the d2 mode and is responsible 
for making the violation of the CS inequality possible. 
Thus, we have proven that the quartic and the quadratic 
violations are equivalent for Id2 , normal-anomalous cor¬ 
relations. 

Now we turn our attention to the GPH criterion. From 
Eq. (1A18D we compute the GPH function for the pair of 
modes Id2: 

'Pld2 = ( gil9d2d2 — |c/d2| 2 )[(<?// + I)(<7d2d2 + 1) ~ | c Jd2| 2 ] • 

( 21 ) 

As the second (square) bracket in the r.h.s. of Eq. 
m is always positive [see Eq. mi we conclude that 
'Pid 2 < 0 if and only if Eq. (usd is satisfied. As the state 
is Gaussian, the GPH criterion is equivalent to the entan¬ 
glement of the state of the system p. We conclude that 
the quadratic and quartic CS violations are equivalent 
to entanglement. We note that the equivalence between 
condition (I19|) and entanglement was already pointed out 
in Ref. [22| but the connection with the CS violation 
criterion of Ref. E2 was not made explicit. A similar 
result appeared in Ref. [45j on the equivalence between 
non-separability and CS violation when studying the cor¬ 
relation function at two different times in fluids of light. 

In regard to the correlation between the two normal 
modes, we obtain, following similar arguments, that the 
quartic CS inequality © reduces to: 

A guugdidi j (22) 

which can never be violated, as previously explained [see 
Eq. (fl2)l and accompanying discussion]. In particular, 
the inequality (E^l) can be rewritten as: 

W u ■ otdi\ 2 < \a u \ 2 \a d i\ 2 , 


(23) 
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which is always satisfied for two complex vectors. Thus, 
there is no CS violation in the correlation between normal 
modes. 

When considering the entanglement between the two 
normal modes, we have \g u di\ > \c u di\ = 0 and thus 
there is no entanglement [see Eq. (1A18I) and the ensuing 
discussion]. 

From the previous discussion we conclude that, when 
considering Gaussian and incoming incoherent states 
[i.e., states satisfying Eq. (ED]. all the here considered 
criteria for characterizing the quantum character of the 
system become equivalent. This result is important be¬ 
cause it unifies the work of Refs. [13, [HI [H[ • 

A most important particular case is that involving a 
thermal distribution of incoming quasiparticles that have 
thermalized in the comoving reference frame [22l . l46j so 
that their occupation factor is: 


exp(fb(w)/T) - 1 


(24) 


with f2j(w) the commoving frequency of the mode i-in 
at laboratory frequency w, as given by Eq. 0. This 
state satisfies the mentioned conditions, so there is no 
difference between using the GPH criterion or the CS 
violation to characterize the quantum behavior of the 
system. 

Finally, we discuss the differences that appear when re¬ 
moving some restrictions. If the state is Gaussian but not 
incoherent in the incoming channels, the GPH criterion is 
still equivalent to the non-separability of the system. On 
the other hand, in the same case (Gaussian and not in¬ 
coherent), the quadratic CS violation is no longer equiv¬ 
alent to the GPH criterion; rather, it is only a sufficient 
condition for it. As a consequence, the GPH criterion is a 
more powerful criterion than the quadratic CS violation 
to detect nonseparability (it can identify nonseparability 
in cases where the quadratic CS violation would fail). 

However, when relaxing the requirement of incoher¬ 
ence, the quartic CS violation becomes independent of 
the GPH criterion (see Appendix lAl). In this context, we 
wish to note, using quantum optics terminology, that the 
presence of quartic CS violation or entanglement requires 
that the system is described by a non-classical Glauber- 
Sudarshan function, i.e., a Glauber-Sudarshan function 
that takes negative values. 

On the other hand, if the state of the system is not 
Gaussian but is incoherent in the incoming channels, the 
quartic CS violation is also independent of the entangle¬ 
ment. The quadratic CS violation Eq. (flTH) and the GPH 
criterion are still equivalent between them, but no longer 
equivalent to the entanglement of the system. 

For a general state which does not satisfy any of the 
previous conditions, the quadratic CS violation is only a 
sufficient condition for the GPH criterion, which in turn 
is a sufficient condition for the presence of entanglement. 
We summarize all these logical relations in Table H) 

Remarkably, the quadratic CS violation Eq. ED re¬ 
veals at the same time two different aspects of quantum 


behavior: the violation of a classical inequality and the 
entanglement of the system. 


IV. ON THE EXPERIMENTAL DETECTION 

In this section, we analyze possible detection schemes 
of the criteria discussed in the previous section. A partic¬ 
ular type of CS violation between two colliding conden¬ 
sates was measured using time-of-flight detection ]34| . 
The possible detection of the quartic CS violation in a 
TOF experiment for analog HR was discussed in Ref. 
E3- A recent work has studied the measurement of quar¬ 
tic CS violations H3| using phonon evaporation, follow¬ 
ing a related work on the dynamical Casimir effect [48j] ■ 
On the other hand, the detection of entanglement us¬ 
ing the GPH criterion was analyzed in Ref. [23] using 
density fluctuations or the optomechanical detection of 
phonons. Recently, it has been proposed that the entan¬ 
glement of Hawking radiation in a BEC can be detected 
experimentally by measuring the density-density correla¬ 
tions through in situ imaging [24j . All these setups have 
in common that they involve collective atom flow and 
thus can be viewed as elementary components of a larger 
atomtronic circuit. 

For illustrative purposes, we focus on the comparison 
between the detection of the CS violation and the GPH 
criterion for the particular scenario of density-density 
correlations, but the analysis which we present in this 
work can be generalized to other schemes. The role of the 
spatial density correlation function in analog models has 
been extensively studied in several works [3, UM [HJ • In 
the BEC context, the measurement of the spatial density 
correlations has been experimentally used to character¬ 
ized the BH Laser [l]|- The density-density correlation 
function also plays an important role in the polariton 
analog [Tol. . 

We start by considering the density in the same BdG 
approximation of Sec. Inland expand up to first order in 
the fluctuations of the field operator: 

n(x) — n 0 (x) + \Jno{x)(j)(x), cj)(x) = ip(x) + fl(x) (25) 

The expression of <j>(x) in terms of the BdG scattering 
states is formally similar to that of tp(x), Eq. 10: 

/•OO _ 

<i>(x) = / du y j rj_i n>a ,(g) 7 j- in (a;) (26) 

Jo j 

/' w max 

+ / d(jJ r d2 -in, W (^)7d2-in( a; ) + Hc * • 

Jo 

with r a ,uj{x) = u a ,w(x) + v a ,w(x). Importantly, we note 
that the density in the BdG approximation is linear in the 
destruction operators, rather than quadratic. Thus, we 
can extract the quadratic correlations of Sec. Ell by mea- 
suring density-density correlations [23, [13 ■ For that pur¬ 
pose, we restrict ourselves to asymptotic (subsonic and 
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supersonic) regions, where the GP wave function adopts 
the form of Eq. ©• We refer to the region between the 
two asymptotic regions as the scattering region, where 
the sonic black hole is placed; i.e., within the scattering 
region the flow velocity crosses the sound velocity at least 
once. In the following, we assume that the scattering re¬ 
gion is placed near x = 0 and its size is much smaller 
than the size of the asymptotic, homogeneous regions. 

By taking the Fourier transform of the density in the 
asymptotic regions at k ^ 0, we can get rid of the con¬ 
densate signal and extract the phonon signal. For def¬ 
initeness, we focus on the correlations u — d2 but the 
procedure for the other cases is similar. We define: 

h r (k) = j d x h{x) f* {x) e~ lkx , (27) 

where f r (x) are normalized functions (f dx\f r (*)| 2 = 1) 
localized in the corresponding asymptotic homogeneous 
regions, sufficiently far from the scattering region. They 
represent the envelope of the Fourier transform, which 


has to be introduced to explicitly take into account the 
fact that the asymptotic subsonic and supersonic regions 
are placed in different spatial regions in a realistic situa¬ 
tion. We choose them such that 

= < 28 > 

with / a symmetric and real dimensionless function, x r 
the point where the envelope is centered, and L r the size 
of the wave packet, which is taken sufficiently large for 
the Fourier transform of the envelope function, 

fr(k) = -^= J dx e~ lkx f r (x) , (29) 

to be sufficiently well peaked at zero momentum. Since 
the subsonic (supersonic) region is placed at the left 
(right) of the BH, we have x u < 0 {xd > 0). 

Taking into account all the previous considerations, we 
obtain: 


n u {k u -ont{u})) ~ f dw' f*(k u -out(u') - fcu-outM)- out^ ) 1/2 7u-out(^') 


|^n—out )| 

+ f*{-k u - in (uj') - k u - outM) . ru ~ m ^ j / 7^-in(k/) 

|w u _in (u/) 

/“tJrnax j, /, ,f\ 

fld{-k d 2-out M) — V™d du/ /d(/Cd 2 -out(w) - fcd 2 -out(w'))-j- , 

Jo \Wd2-out (<*/) | 1/2 


7d2-out(w') 


+ J do;' fd(kd2-out{u) - k dl - out (uj')) rdl ou / t ^| 1/2 7di- ou t( u,, )» 


|Wdl-out (<*/) | 1/2 


(30) 


where r a (w) is [see Eq. ©] 

r a (w) = U a (u>) + v a {io ) = 


k 2 n 


(31) 


2 fi Q (w)’ 

and we have used \u> — v r k a {oj)\ = $7 a (u;) = \Jc^k 2 + 


Now, we can connect the correlations studied in Sec. m 
with the density correlations by taking into account that 
h\(k) = n r (—k). We obtain: 


— {flu ( k u — ou t (^)) f'U {k U — out (^))) — Xl u T u _ out (^) \9uu{l*f) T ^w(^uu(w)) “t“ 1] 
G d 2d2(w) = {fid (fc d 2-out(w)) fl d {~k d 2-out M)) = ^d^-out ( W ) [9d2d2 M + 9 dldl{^ dld2^)) + 1] 
G ud2{^') = {fid ( kd2— out (^0) H u {ku— out (^0)) — F(cu)-y/ Tl u TldTu— out l'd2—out^ud2 (^) ) 


(32) 


I- 

with Wdid 2 (w) = Wdi(fcd 2 -out(w)) and co uu ( w) = function between the subsonic and supersonic regions, 
Uu{—k u - out(w)), Wj(fc) being the dispersion relation of 
the mode i. Their corresponding values are obtained 
graphically in Fig. [2] On the other hand, the overlap 
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F(ui), is given by: 


F(u) = / dfc /: 

i r 

dxf* 


y/ L u Ld 


C(w 

£ (u>) X + X. 


fd (*CM) 


(33) 


/* 


C M 1 z - 

Ld 


with 


CM 


yju —out (cu) 

Wd2-out(w) 


1/2 


(34) 


The previous integral can be interpreted as the scalar 
product of two normalized functions, which satisfies 
F(u) <1. It is easy to prove that F(u) = 1 when 


CM = 

This condition has been interpreted as the condition for 
the wave-packets to have equal width in frequency space 
E2- However, resulting from the well-known properties 
of the spatial density correlations, Eq. © can be also 
interpreted as expressing that the envelopes in the sub¬ 
sonic and supersonic regions have to be placed along the 
correlation lines that maximize the spatial density corre¬ 
lation function, see Refs. 01110 - 

In the following, we suppose that condition (l35l) is sat¬ 
isfied and F(uj) = 1. It is worth noting that the correla¬ 
tion functions Gij{uj) can never violate the CS inequality 
by themselves since G u d 2 < \ZG U uGd 2 d 2 ■ Thus, we try 
to relate these functions, which can be measured, to the 
phonon correlations considered in the previous section. 
For that purpose, we normalize the density-density cor¬ 
relations of Eq. (l32l) : 



Gij ( w ) 


_MM_ 

yjTliTljTi —out (a^)rj'—out 


(36) 


where n t ,j must be interpreted as = rid when the 
indices i. j take values dl,d2. From Eq. (13^1) . we see 
that Gud 2 {w) = c u d 2 (u). For the extraction of the other 
correlations, in a similar way to Ref. (24j, we define the 
magnitudes: 


9uu{^') — ^liu(^) T ^u(^u«(^)) I 

£?<22<22M — 9d2d2^) + Sdldl Mld2M) = £<22(22 M ~ 1 

(37) 


The —1 appearing in the definitions of g U u{w), gd 2 d 2 (w) 
reflects the subtraction of the atomic depletion contribu¬ 
tion. 

We see that g U u,gd 2 d 2 are over-estimations of the cor¬ 
relation functions g uu ,gd 2 d 2 , since n u ,gdidi are always 
positive. From this, we can define 


A u d 2 =Jc u d 2 \ 2 — g uu gd 2 d 2 , _ (38) 

noting that A„ d2 > A u d 2 - Thus, measuring A, ld2 > 0 
amounts to experimentally observing the quadratic CS 
violation. We note that, for these calculations, the sole 
assumption has been made that the state of the system 
is stationary. 

The upshot of this discussion is that we can observe 
the violation of a quadratic CS inequality through the 
measurement of the function A u d 2 - Similar claims can 
be made about the measurement of the quadratic CS 
inequality involving the dl, d2 modes. We can repeat 
the same strategy as before to obtain the correlation 
functions gdidi,Cdid 2 - However, it is not possible to ob¬ 
tain the correlation functions c uu ,c d idi,Cd 2 d 2 ,g u d 2 ,gdid 2 
from this scheme. This is due to the vanishing over¬ 
lap integrals that appear when trying to obtain the cor¬ 
responding correlations. For instance, if we compute 
{hu {k u - o U t (u))n u {k u - out M)), we face an overlap in¬ 
tegral of the type 


I d k /; (fc) /: (fc) = Jdx f: (x) f* (-*) = 0, (39) 


because the f u function is, by construction, well local¬ 
ized in the subsonic region, far from the scattering region 
around x = 0, see Eq. ED and the discussion below. 
Thus, we cannot obtain c uu by this method. A similar 
reasoning applies to other correlations. 

It is important to remark on the crucial role played by 
the spatial location of the asymptotic regions. If we did 
not introduce the envelopes in Eq. ED and rather did 
take the Fourier transforms as ideally infinite, we would 
obtain: 


J 


Hu {ku —out (^)) — \/ 

ftd(-fcd2-outM) - V™d 


T u —out (m) 

|w u -out Ml 1/2 


Tu—out(^) H - 


T it—in MmM) 
^u —in 0^)) 


7;_mM«M) 


f ^ d2— out(^) ^ / \ . ^ dl— out {^dld2 (^)) ~ f / / xx\ 

[\wd 2 -out ( W )|i/2 7d2 -° ut(w) + Mr-out M ld2 M) |l /2 7«il-outMl«i2(a;))J , 


( 40 ) 























k k 

FIG. 2: Graphical obtention of L 0 did 2 (u) and —oj uu (lo). See Fig. [T]for mode notation. 


We can infer that such an approach would lead to contradictory results. For instance, the commutator 

[n u (k u - ou t(uj)) , h u (-fc u _ in (w))] = n u — 7 “•-° llt ( a; ) — ? u-mM s uu (w) ± 0 (41) 

K-out HI 7 K-in H | 1/2 

would be non-vanishing, in contradiction with the well-known fact that it has to be zero since [h(x),n(x')\ = 0 for 
every x, x'. By way of contrast, when using the full expressions of Eq. (I30|) we arrive at the correct result 


\flu (^ii—out(k-O) 5 Tlu ( ku— in(^))] — Tin 


T u —out (^) 


i— in(^) 


H- ou t Hl 1/2 Hl 1/2 

i- 


s„ 


t (w) J do/ f* 


\w u — in (k-0 


/: - 


\lX U —out (&) 


i42) 


The integral of the r.h.s, after returning to real space, is 
similar to that considered in Eq. (15^)) and gives zero. 

We arrive at the same conclusion when consider¬ 
ing the density-density correlations in real space. As 
noted in Ref. [20j], only the terms with a station¬ 
ary phase should be kept when computing (n(x)fi(x ')), 
which yields the condition (1351) . In particular, only 
the “proper” correlations between out-out modes can 
be extracted c u d 2 , Cdid 2 , gudi- The other correlations 
Cuu,Cdidi,c d2d 2,gud2,gdid2, c u di , even when they are non¬ 
zero, cannot be obtained because the associated expo¬ 
nential terms do not present a stationary phase when 
integrating over frequencies. 

We conclude that an important consequence of tak¬ 
ing into account the spatial location of the subsonic and 
supersonic regions is that only a limited number of cor¬ 
relations can be experimentally observed. In particular, 
this implies that we cannot measure all the correlation 
functions appearing in the GPH criterion (IA11I) by this 
procedure. Moreover, and for the same reason, if the 
state is Gaussian, we cannot obtain all the correlation 
functions appearing in the quartic CS violation, see Eq. 

EnD. 

The results of the present section show that, within 


this kind of detection schemes, we can only aim at ob¬ 
serving a quadratic CS violation. One can expect this 
limitation to be not exclusive of the specific procedure 
here considered. The reason is that any realistic attempt 
to obtain the correlation functions between phonons must 
necessarily take into account the spatial location of the 
asymptotic regions, which implies that a similar reason¬ 
ing does apply. 


Nevertheless, the previous considerations do not pose a 
problem for the detection of the quantum Hawking effect 
for two reasons: (a) if the state of the system is Gaus¬ 
sian and incoherent over the incoming channels, we have 
proven that the GPH criterion and the quartic CS vio¬ 
lation are equivalent to the quadratic CS violation and 
(b) even if the state of the system does not belong to 
that class, the quadratic CS violation is still a signature 
of the presence of the entanglement; in particular, it is 
a sufficient condition for the fulfillment of the GPH cri¬ 
terion. Even more than that, as remarked at the end of 
Sec. EH the quadratic CS violation is also by itself a 
clear indication of the quantum nature of the system, as 
CS inequalities are always satisfied in a classical system. 
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V. NUMERICAL RESULTS 

We investigate here the possibility of experimental de¬ 
tection of CS violation in different black-hole setups. For 
that purpose, we compute the measurable quantity A ud 2 , 
defined in Eq. (1551) . We suppose that the state of the sys¬ 
tem is given by a thermal distribution in the incoming 
modes, as expressed in Eq. (l24l) . Since we only focus on 
the identification of some physical trends and not on the 
study of the whole parameter space, characterized by 7 
variables (see Appendix [Bj, we will consider for simplic¬ 
ity some typical configurations studied in the literature in 
order to compute A ud 2 - For this analysis we distinguish 
between non-resonant and resonant structures. For the 
non-resonant case, we consider in this work two scenarios: 
the single delta barrier configuration and the waterfall 
configuration, schematically depicted in Fig. [3] 

In the delta barrier configuration, the black hole forms 
near a localized potential of the form V(x) = Z5(x) 
mum-, see left panel of Fig. [3] This configuration 
permits us to study theoretically the flow of a condensate 
through a potential barrier of finite size. On the other 
hand, the waterfall configuration creates the sonic hori¬ 
zon by using a negative step potential V(x) = — Vq6(x), 
where 9(x) is the step function . This kind of scenario, 
where the black hole is produced by introducing a nega¬ 
tive potential, has been already experimentally realized 
by the Technion group [12]. 

In Fig. 01 we represent the Hawking spectrum, \S u d 2 \ 2 
and the functions A u d 2 , A u d 2 for the delta barrier and 
the waterfall setups, both at temperature T = 0.2, which 
is of the order of magnitude of typical experimental se¬ 
tups, where it can be as low as T ~ 0.1 [ll]. We see 
that, for some range of frequencies, the measurable func¬ 
tion A u d 2 satisfies A u d 2 > 0, which implies the presence 
of CS violation. As noted previously, A ud 2 > 0 is a more 
restrictive condition than the bare quadratic CS viola¬ 
tion A u d 2 > 0. However, we observe that for sufficiently 
large w, both A ud 2 ,A ud2 converge to 15^21 2 , since for 
high ui we have that the occupation numbers are negligi¬ 
ble and thus, we recover the zero-temperature limit, for 
which A ud 2 (T = 0) = IS^I 2 - Moreover, we see that, at 
high uj, \Sdid 2 {u d id 2 {w))\ 2 = 0, because for ui such that 
Wdid 2 (w) > w max the anomalous scattering channel disap¬ 
pears, and then we have A ud2 ~ A ud2 (T = 0) = |S„d 2 | 2 - 
Finally, we note that there is no violation near w = 0, as 
argued quite generally in Ref. [17]. 

We now focus on resonant configurations. In partic¬ 
ular, we focus on the case of a double delta barrier, 
which was first considered in Ref. [l3] . In this situa¬ 
tion, the black hole is formed by two single delta barriers 
separated by a distance d, so the potential is given by 
V(x) = Z[5(x) + S(x — d)}. This setup is schematically 
depicted in Fig. [5] It was shown in Ref. 0 that reso¬ 
nant spectra can be expected to present a strong signal of 
CS violation. This trend can be observed in Fig. [Gl where 
we represent the same magnitudes as in Fig. 0]but now 
for a higher temperature T = 0.7. We see that, even for 


this relatively high temperature, the experimental signal 
A ud 2 is, at the resonance frequency, substantially larger 
than in the non-resonant case. 


VI. CONCLUSIONS 

We have analyzed the range of validity of different ex¬ 
istent criteria for the identification of deep quantum be¬ 
havior as applied to the outgoing channels involved in 
the spontaneous Hawking radiation. Specifically, we have 
compared the generalized Peres-Horodecki criterion with 
the violation of quadratic and quartic Cauchy-Schwarz 
inequalities. We have shown that, under certain physi¬ 
cal conditions (Gaussian processes and, simultaneously, 
incoherent incoming channels), all the considered crite¬ 
ria are equivalent. When such conditions are not ful¬ 
filled, we have shown that the quartic CS violation and 
the non-separability of the state represent independent 
mathematical conditions. 

We have also investigated the possible measurement 
of the different criteria in realistic scenarios. By taking 
into account the different spatial location of the subsonic 
and supersonic regions, we have shown that only certain 
correlation functions can be obtained. For simplicity we 
have focused only on one detection scheme, but we expect 
similar problems to arise in other kind of measurements, 
since in any realistic situation the supersonic and sub¬ 
sonic regions are necessarily placed in different regions. 
However, our work also shows that this limitation is not a 
major problem, as we can often measure the quadratic CS 
violation, which is also a sufficient condition for the en¬ 
tanglement of the system. Finally, our numerical results 
of Section IVl show that, in typical analog configurations, 
the CS violation can be detected in an achievable range 
of temperatures. 

Sonic black holes, as well as related setups such as 
black-hole lasers, have the potential to become an im¬ 
portant element of atomtronic circuits, for instance, as a 
source of entangled quasiparticles in the broader context 
of quantum communication. 
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FIG. 3: Scheme of the delta barrier configuration (left) and the waterfall configuration (right). Here, c(x),v(x) are the local 
sound and flow velocities, while V (®) is the external potential. 




FIG. 4: Plot of the Hawking radiation spectrum and the functions A u d 2 , A„d 2 for a delta barrier configuration (left panel) 
and a waterfall configuration (right panel), at temperature T = 0.2. For the delta barrier, the subsonic flow speed is v u = 0.3 
and the barrier strength is Z = 0.62. For the waterfall setup, the subsonic flow speed is v u = 0.5 and the potential depth is 
Vo = 1.125. The corresponding Hawking temperatures [22l. [46] are Th = 0.24 and Th = 0.14, respectively 


Appendix A: General remarks on Cauchy-Schwarz 
inequalities and entanglement 

1. Cauchy-Schwarz inequality for operators 


We briefly review the mathematical CS inequalities for 
operators in a Hilbert space, which are always fulfilled, 
see Ref. [43j for example. We start by defining the scalar 
product associated to a state p for two operator A, B as: 


(A, B) = (A*B) = T 'r(pA'B). (Al) 

It is easy to check that the previous product satisfies 
the usual properties of a scalar product. Thus, it satisfies 
the mathematical Cauchy-Schwarz inequality: 

| (A^B) I 2 < (i f i) (&B) . (A2) 

In particular, the previous inequality implies for B = 1: 

|(i)| 2 = |(it)|2< ( iti ) , (A3) 
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which also leaves the possibility of quartic CS violations; 
see Eq. (fTUl) . 


2. Quadratic Cauchy-Schwarz violation and the 
generalized Peres-Horodecki criterion 

We now consider bipartite states of two modes i.j. We 
define the vector X = <7j,j5j] T , where the phase 

space operators qk,Pk are related to the annihilation op¬ 
erator of the modes k = i,j through 7 /. = (g*, +ipk)/V 2 . 
Following Ref. [42j], we define the 4x4 matrix: 


FIG. 5: Scheme of the double delta configuration. 



0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 

& / ^max 


FIG. 6: Same as Fig. [4] but now for a double delta bar¬ 
rier configuration with parameters Z = 2.2, d = 3.62 and 
v u = 0.01. The temperature of the system is T = 0.7. The 
inset zooms into the peak region. In this case, the Hawking 
temperature of the two horizons (see Fig. is T h « 0.13. 


and, for A = 7 * and B = 7 j, 

\9ij\ 2 = I (7i7j> I 2 < (7?7») (tfjli) = 9a9jj • (A4) 

Thus, the CS inequality \gij\ 2 < gugjj is always satisfied 
and can never be violated. On the other hand, taking 
A = 7 J and B = Jj, we arrive at: 

|cij| 2 = | {lilj) I 2 < {lill) (lj7i) = (9u + 1 )9jj • (A5) 

Thus, the mathematical CS inequality states \cij\ 2 < 
(gu + 1 )gjj, leaving the possibility of quadratic CS vi¬ 
olation, Eq. m- The same type of argument holds 
for quartic CS violations, where an analogous reasoning 
leads to: 

l r ij| 2 = I (7i7i7]7j) I 2 < (lllillli) iljljljlj) 

= (Tu + gu){Tjj + gjj) (A6) 


M = V + i— (A7) 


w = 

J 

0 

, J = 

0 1 


0 

J 


-1 0 


where V is the covariance matrix with matrix elements 
V a p = \{XX a , XXp}, with AX = X — (X) and 
{AX a , AXg} the anticommutator. The matrix J is the 
simplectic matrix in two dimensions. The matrix M is 
non-negative, M > 0, as this fact represents an alterna¬ 
tive expression of the uncertainty principle [42| . 

If the state p is separable, it has the form (fl3l) . and 
by taking partial transpose of the state with respect to 
the subsystem j, we obtain an operator p t which is also 
necessarily a physical state. The GPH criterion is based 
on this fact. It can be shown that if pt is a physical state, 
the matrix M t is also non-negative, where M t is given by 
0 : 


M t = V t + i— (A 8 ) 

V t = AEA, A = diag[l, 1,1, —1]. 

The condition M t > 0 is the uncertainty principle for the 
state pt- As M is non-negative, M t is non-negative if 
and only if det M t > 0. The conditions det M, M t > 0 
are equivalent to V± > 0 , respectively, where: 


= det A t det Aj + ^ — 7 det CJ r j (A9) 

— tr (JAiJCijJAjJCjj) — i(det A, + det Aj). 

The matrices Aj, Aj , Cjj are 2x2 submatrices of the 
covariance matrix V : 


Aj Cjj 

^5 


(A10) 


We note that V+ > 0 is always satisfied since p is a 
physical state. However, when V~j < 0, the state p t is 
not a physical state, which implies that the original state 
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p is not separable. We can put together both conditions 
by defining the GPH function Vij as: 

Vij = det Ai det Aj + ^ — | detCjj|^ (All) 

— tr (J AiJCijJ AjJCjj) — i(det Ai + det Aj), 

Thus, if Vij < 0, the state is entangled. This result 
is the GPH criterion. However, when det Cij > 0, the 
state is separable, as Vij = > 0, so only states with 

det C\j < 0 can be entangled. 

On the basis of the previous results, we now prove 
that the quadratic CS violation implies the fulfillment 
of the GPH criterion. Suppose that Vij > 0. In that 
case, the matrix M t is non-negative and we can define an 
associated scalar product for vectors u, v £ C 4 as: 


By inserting 



1 

O 

1_ 


’ 1 " 

1 

0 

1 

i 

u = —j= 

V2 

1 


0 


i 


1 

O 

_1 


into Eq. (1A13I) . we obtain the quadratic CS inequality 
\c-ij \ < gugjj. Thus, if there is quadratic CS violation, 
the matrix M t cannot be non-negative, which implies 
that Vij < 0 and the GPH criterion is satisfied. We 
conclude that the quadratic CS violation is a sufficient 
condition for the fulfillment of the GPH criterion. 


3. Quartic Cauchy-Schwarz violation and 
entanglement 


(u,v)t = u*M t v (A12) 

which satisfies the associated CS inequality: 

\(u,v) t \ 2 <(u,u) t (v,v) t . (A13) 

Assuming the usual case where the annihilation and de¬ 
struction operators have zero expectation values, the ma¬ 
trices Aj j and Cij can be written in terms of the corre¬ 
lation functions of Eq. m as: 


Ak = 


Cij = 


9kk 


I 2 


He Cfzk Im c 
Im Cfefc -Re c k k 


He (fjij Cjj ) Ini (fjj,j -f- Cij ) 

Im {-gij + c^) R e(gij - dj) 


(A14) 

(A15) 


with the index k = i,j and I 2 the 2x2 identity matrix. 


The previous arguments cannot be applied to the quar¬ 
tic CS violation Eq. (in. As a counterexample, the di¬ 
rect product of two pure number states of the modes i , j. 
p = \n) (n| ® |n , )(ro , |, which is a manifestly separable 
state, violates the quartic CS inequality. Even for Gaus¬ 
sian states, the two conditions, quartic CS violation and 
entanglement, are still independent. For instance, for a 
general Gaussian state, the quartic CS violation of Eq. 
ms can be expressed, via Wick’s theorem, in terms of 
the quadratic correlations: 


— \ c ij\ T gugjj tjtyii + \ C ii\ 2 ^^djj + \ c jj \ 2 

(A17) 

For Gaussian states the GPH condition is equivalent 
to the non-separability of the system. Using Eqs. (IA11I) . 
(IA14I) . we can compute explicitly the GPH function: 


J 


'Pij — [dagjj \ c ij\ }[{gu + C)(gjj + 1) \cij | ] + 4 I gu + I R e (gij c jjCij) + 4 I gjj + ) Re(gijCac i A 


1 


2 j ■ - \ jjj ■ 2 

2 1 |„ |4 \„ 12 / „ 1 „ 1 „ 1 o I |2\ 


1 


2Re(c i jC ii Cjj) 2Re(g i jCnCjj) + \cu\ \cjj\ +|gij| \gij\ (gu + 9jj + 2gagjj + 21c*^| ) 
~ gu(gu + l)l c jj| ^ + l)l c ii| • 


(A18) 


This expression applies whenever \cij\ > \gij\- For \cij\ < 
| g^ |, we have det Cij > 0 and the state is separable, as 
previously explained. 

As we can see, the expressions for the quartic CS viola¬ 
tion and the GPH criterion represent different conditions 
for arbitrary Gaussian states. In particular, it is easy 
to find states violating the quartic CS inequality with 
det Cij > 0, which means that they are separable. Also, 
we can find entangled states that satisfy the quartic CS 
inequality. 


Appendix B: Parametrization of outgoing 
correlations for incoherent incoming Gaussian states 


We discuss in this Appendix the parameters needed 
to compute the correlation functions defined in the main 
text. A related discussion appeared in Ref. [22j. The 
scattering matrix S relates the “out” scattering states 
with the “in” scattering states through Eq. (Q. The 
conservation of the commutation relations implies the re- 
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lation: 


S^rjS = r] = diag(l, 1,-1), (Bl) 


which means that S 6 U( 2,1). It follows that S 1 = 
rjS'r) and, using standard linear algebra, we can prove: 


n* _ 111 1 i 'll. I 

ij ~ det S 


(B2) 


where rriij is the minor associated to the S matrix ele¬ 
ment Sij. In order to simplify the notation, in this section 
we relabel the indices it, dl, d2 as 1, 2,3 in order to match 
the matrix indices ordering. With the help of the previ¬ 
ous results, it can be proven that any matrix S € 1/(2,1) 
satisfies: 



det S 

0 

0 

BSA = Spa = 

0 

c* 

°33 

N 


0 

N 

<§33 

1 

-1 

1 

Co 

to 

CO 

Sl 3 

0 

b =n 

Q* 

°13 

c* 

°23 

0 

0 

0 

N 


-S32 <§ 3 \ o 
<§32 0 
0 0 JV 


(B3) 


where N = v / |S'i 3 | 2 + |<§ 2 3 | 2 = = 

x/l'S'as | 2 — 1- We see that the matrices A, B are unitary. 
Inverting the previous relation, we obtain: 


S = B^S pa AK (B4) 


Equations mm show explicitly that we only 
need 9 real parameters to parametrize the S matrix: 


|<§3l|j |<§1 3 |; |<§3 3 |)Xi 013i 023j‘/’SI, 032, 033) W ^h fiij the 
phase of Sij, Sij = \Sij\e z ^ i: > and x that of det 5, 
det 5 = e lx . That is, we use 3 amplitudes and 6 phases 
to characterize the complete scattering matrix. Interest¬ 
ingly, the matrix Spa shows clearly that the Hawking 
radiation acts as a non-degenerate parametric amplifier 
& 

A convenient way to write the S matrix for our problem 
is: 

s = Ul t SU in (B5) 

U k = diag[e i5l_fc , e iS2 ~ k , e i&3 ~ k ] 

(with k = in, out), since we can fix the elements of 
U in . ou t in such a way that they absorb the phases 

032, ^13, 023,033, and so <§ 3 i, § 32 , <§ 13 , § 23 j <§33 be¬ 
come purely real. This means that the matrix S is char¬ 
acterized by only 4 parameters, |<§ 3 i|, |<S'i 3 |, |<S 33 |,x, with 

det § = e ix = det S e - i ^ 13+ ^ 3+ ^ 13 +^ 3 -^ 3 ). (B6) 


Nevertheless, we do not need all 9 parameters for com¬ 
puting the physical quantities of the main text. For in¬ 
stance, the CS violations of Eq. m are invariant under 
phase transformations of the “out” states. It can be also 
shown that the GPH function is invariant under these 
transformations, see Eq. (IA18I) . As we are always con¬ 
sidering incoherent incoming modes, all the quantities are 
also invariant under phase transformations of the incom¬ 
ing modes. From Eq. (IB5I) . and taking into account the 
previous observations, we conclude that all the requested 
quantities depend only on the matrix S , which is char¬ 
acterized by just 4 parameters, as noted in the previous 
paragraph. Using Eqs. (IBljl . (1B4I) we can write § as: 



— sin a 

cos a 

0 


e ix 

0 

0 


— sin/? 

cos /? 

0 


S = 

COS CM 

sin a 

0 


0 

coshy 

sinh y 


cos /? 

sin/3 

0 

(B7) 


0 

0 

1 


0 

sinh y 

coshy 


0 

0 

1 














with |<S 33 | = coshy, |Si 3 | = sinhycosa, and |<S 3 i| = 
sinhycos/3. Thus, we only need 4 parameters in order 
to parametrize the scattering matrix for the calculations 
of this work: a , /3, y, We note that this parametriza- 
tion differs slightly from that used in Ref. [HJ, where 
4 amplitudes were used to characterize the state of the 
system. Nevertheless, that election is equivalent to that 


used here. 

Finally, taking into account that the state of the sys¬ 
tem p is characterized by only 3 numbers, n, [with 
i = l,2,3, see Eq. m\, we have that the whole problem 
is completely determined by 7 parameters, 4 arising from 
the S matrix and 3 arising from the specification of the 
incoherent incoming Gaussian state. 
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